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Beschreibung 

Verfahren zur Grauwert-basierten Bildf ilterung in der Compu- 
ter-Tomographie 

Mit modernen medizinischen Diagnoseverf ahren, wie beispiels- 
weise der Rontgen-Computertomographie (CT) , konnen Bilddaten 
eines untersuchten Messobjektes gewonnen werden. In der Regel 
handelt es sich bei dem untersuchten Messobjekt um einen Pa- 
tienten . 

Die Rontgen-Computertomographie - im nachf olgenden kurz mit 
CT bezeichnet - ist ein spezielles Rontgen-Auf nahmeverf ahren, 
dass sich im Bildaufbau grundsatzlich von dem klassischen 
Rontgen-Schichtauf nahmeverf ahren unterscheidet . Bei CT- 
Aufnahmen erhalt man Transversal-Schnittbilder , also Abbil- 
dungen von Korperschichten, die im wesent lichen senkrecht zur 
Korperachse orientiert sind. Die im Bild dargestellte gewebe- 
spezifische physikalische Gr6£e ist die Verteilung des Schwa- 
chungswertes von Rontgenstrahlung fi(x,y) in der Schnittebene . 
Das CT-Bild erhalt man durch Rekonstruktion der vom verwende- 
ten MeSsystem gelieferten eindimensionalen Projektionen der 
zweidimensionalen Verteilung von |a(x,y) aus zahlreichen ver- 
schiedenen Blickwinkeln . 

Die Projektionsdaten ermitteln sich aus der Intensitat I ei- 
nes Rontgenstrahls nach seinem Weg durch die abzubildende 
Schicht und seiner ursprunglichen Intensitat I 0 an der Ront- 
genquelle gemafi dem Absorptionsgesetz 



Der Integrationsweg L reprasentiert die Bahn des betrachteten 
Rontgenstrahls durch die zweidimensionale Schwachungsvertei- 
lung |X(x,y) . Eine Bildpro j ektion setzt sich dann aus den mit 




(1) 
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den Rontgenstrahlen einer Blickrichtung gewonnenen Messwerten 
der Linienintegrale durch die Objektschicht zusammen. 

Man erhalt die aus verschiedensten Richtungen stammenden Pro- 
5 jektionen - charakterisiert durch den Pro jektionswinkel a - 
durch ein kombiniertes Rontgenrohren-Detektor-System, das in 
der Schichtebene um das Objekt rotiert. Die derzeit gebrauch- 
lichsten Gerate sind sogenannte „Facherstrahlgerate" bei de- 
nen Rohre und ein Array aus Detektoren (eine lineare Anord- 

10 nung von Detektoren) in der Schichtebene gemeinsam um ein 

Drehzentrum, welches auch Mitte des kreisf ormigen Messfeldes 
ist, rotieren. Die mit sehr langen Messzeiten behafteten „Pa- 
rallelstrahlgerate" werden hier nicht erlautert . Es sei je- 
doch darauf hingewiesen, dass eine Transformation von Facher- 

15 auf Parallelpro jektionen und umgekehrt moglich ist, so dass 

die vorliegende Erfindung, die anhand eines Facherstrahlgera- 
tes erklart werden soil, ohne Einschrankung auch fur Paral- 
lelstrahlgerate anwendbar ist. 

2 0 Bei Facherstrahlgeometrie besteht eine CT-Aufnahme aus Li- 

nienintegralmesswerten -ln(I/I 0 ) eintref fender Strahlen, die 
durch eine zweidimensionale Verknupfung des Pro jektionswin- 
kels ae[0,2n) und den die Detektorposi tionen def inierenden 
Facherwinkeln Pe[-|3o/Po] (Po ist der halbe Facherof fnungswin- 
^ts kel) charakterisiert sind. Da das MeSsystem nur uber eine 

endliche Anzahl k von Detektorelementen verfiigt und eine Mes- 
sung aus einer endlichen Anzahl y von Pro jektionen besteht 
ist diese Verknupfung diskret und kann durch eine Matrix dar- 
gestellt werden: 

30 P(a y ,&): [0,27c) x [-Po,Po] (2) 

bzw. 



35 



p{y,k): (1,2,... N P ) X (1,2,... N s ) (3) 
Die Matrix p ( y ,k) heifit Sinugramm fiir Facherstrahlgeometrie. 
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Die Pro jektionszahl y und die Kanalzahl k liegen in der Gro- 
Senordnung von 1000. 

Bildet man gemaE Gleichung (1) die Logarithmen, so erhalt man 
5 also die Linienintegrale samtlicher Projektionen 



deren Gesamtheit auch als die Radon-Transf ormierte der Ver- 
10 teilung ^i(x,y) bezeichnet wird. Eine solche Radon- 
L Transforation ist umkehrbar, |4,(x,y) demnach aus p(a,0) durch 
^ Rucktrans formation (inverse Radontrans formation) errechenbar. 
Bei dieser Riicktransf ormation wird ublicherweise ein Fal- 
tungs-Algorithmus angewendet, bei dem die Linienintegrale je 
15 Projektion zunachst mit einer speziellen Funktion gefaltet 
und dann langs der ursprunglichen Strahlrichtungen auf die 
Bildebene ruckproj iziert werden. Diese spezielle Funktion, 
durch die der Faltungs-Algorithmus im wesentlichen charakte- 
risiert ist, wird als "Faltungskern" bezeichnet. Durch die 
20 mathematische Gestalt des Faltungskerns besteht die Moglich- 
keit, die Bildqualitat bei der Rekonstruktion eines CT-Bildes 
aus den CT-Rohdaten gezielt zu beeinf lussen . Beispielsweise 
konnen durch einen entsprechenden Faltungskern hohe Frequen- 
zen betont werden, um die Ortsauf losung im Bild zu erhohen, 
25 oder durch einen entsprechend anders gearteten Faltungskern 

hohe Frequenzen zur Reduktion des Bildrauschens gedampft wer- 
den. Zusammengef asst kann also gesagt werden, dass bei der 
Bildrekonstruktion in der Computertomographie durch Wahl ei- 
nes geeigneten Faltungskernes die Bildcharakteristik, die 
3 0 durch Bildscharf e/Bildkontrast und Bildrauschen (beide ver- 

halten sich zueinander komplementar ) charakterisiert ist, be- 
einf lusst werden kann. 




35 



Auf das Prinzip der Bildrekonstruktion in der CT durch Be- 
rechnung der |i,-Wert-Ver teilung soil nun nicht weiter einge- 
gangen werden. Eine ausfuhrliche Beschreibung der CT- 
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Bildrekonstruktion ist beispielsweise in „Bildgebende Systeme 
fiir die medizinische Diagnostik" , 3 . Auf lage, Miinchen: Publicis 
MCD Verlag, 1995, Hrsg. : Morneburg Heinz, ISBN 3-89578-002-2, 
dargestellt . 

Allerdings ist mit der Berechnung der |i-Wert~Verteilung der 
durchstrahlten Schicht die Aufgabe der Bildrekonstruktion 
noch nicht abgeschlossen . Die Verteilung des Schwachungskoef - 
fizienten ]X reprasentiert im medizinischen Anwendungsbereich 
nur eine anatomische Struktur, welche noch in der Form eines 
Rontgenbildes dargestellt werden muss. 

Nach einem Vorschlag von G.N. Hounsf ield ist es allgemein ub- 
lich geworden, die Werte des linearen Schwachungskoef fizien- 
ten jo, (der die Ma£einheit cm' 1 hat) auf eine dimensionslose 
Skala zu transf ormieren, in der Wasser den Wert 0 und Luft 
den Wert -1000 erhalt. Die Umrechnungsf ormel auf diese „CT- 
Zahl" lautet: 

CT - Zahl = — f^SSSSL 1QQ0 (4) 
/^Vasser 

Die Einheit der CT-Zahl heiSt „Hounsf ield-Unit " (HU) . Diese 
als "Houndsf ield-Skala" bezeichnete Skala ist sehr gut zur 
Darstellung anatomischen Gewebes geeignet, da die Einheit HU 
die Abweichung in Promille von passer ausdriickt und die ^~ 
Werte der meisten korpereigenen Substanzen sich nur wenig vom 
(X-Wert des Wassers unterscheiden. Aus dem Zahlenbereich (von 
-1000 fiir Luft bis ca . 3 000) werden nur ganze Zahlen als Tra- 
ger der Bildinf ormation verwendet . 

Allerdings wurde die Darstellung des gesamten Skalenbereiches 
von etwa 4000 Werten das Unterscheidungsverm6gen des mensch- 
lichen Auges bei weitem iibersteigen. Zudem interessiert den 
Betrachter oft nur ein kleiner Ausschnitt des Schwachungswer- 
tebereiches, z.B. die Dif f erenzierung von grauer und weiSer 
Gehirnsubstanz , die sich nur urn etwa 10 HU unterscheiden. 
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Aus diesem Grunde benutzt man die sogenannte Bildf ensterung 
(engl.: Windowing). Dabei wird nur ein Teil der CT-Werteskala 
ausgewahlt und iiber alle verfiigbaren Graustufen gespreizt. 
Auch kleine Schwachungsunterschiede innerhalb des gewahlten 
5 Fensters werden so zu wahrnehmbaren Grautonunterschieden, 

wahrend alle CT-Werte unterhalb des Fensters schwarz und alle 
CT-Werte oberhalb des Fensters weifi dargestellt sind. Das 
Bildf ens ter kann sowohl in seinem Zentralniveau als auch in 
seiner Weite beliebig variiert werden. 

10 

Nun ist es in der C omputer-Tomo graph ie von Interesse bei be- 
^ stimmten Aufnahmen organspezif ische Einstellungen der Bild- 
charakteristik sowie u.U. ein organspezif isches Windowing 
vorzunehmen. Beispielsweise bei Transversal-Schnittbildern 
15 durch den Brustraum - in welchen gleichzeitig Herz, Lunge, 

Wirbelsaule aufgenommen werden - fiihrt eine organspezif ische 
Optimierung der Bilddarstellung zu einer weitaus besseren U- 
bersicht und erleichtert dem Anwender die Interpretation der 
CT-Auf nahme . 

20 

Eine derart optimierte Aufnahme erfolgt nach dem Stand der 
Technik dadurch, dass nach Aufnahme der betreffenden Schicht 
durch Benutzung unterschiedlicher Faltungskerne bei der Bild- 
rekonstruktion aus den Rohdaten eine Serie von Bildern er- 
zeugt wird die sich jeweils durch unterschiedliche Bildcha- 
rakteristik (Kontrast, Rauschen) voneinander unterscheiden. 
Der Anwender entscheidet dann in welchem Bild das jeweilige 
Organ entsprechend der diagnostischen Anforderung optimal 
dargestellt ist. In den ausgewahlten Bildern muss der Anwen- 

3 0 der das jeweilige Organ segmentieren - in anderen Worten: 

"Markieren" und " Ausschneiden" - und in das endgultige End- 
Bild einfiigen. Zum Segmentieren stehen dem Anwender sogenann- 
te Segmentier-Algorithmen zur Verfugung. Diese f unktionieren 
zumeist derart, dass innerhalb des zu segmentierenden Organs 

3 5 durch den Anwender ein Startpunkt gesetzt wird von dem aus 
nach unterschiedlichen Gesichtspunkten der Rand des Organs 
ermittelt wird. Der Algorithmus hangelt sich an der Organ- 
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grenze entlang bis das ganze Organ abgetastet ist und somit 
ausgeschnitten und in das End-Bild eingefugt werden kann. 

Das Vorgehen bei einer derartigen Segmentierung nach dem 
5 Stand der Technik ist sehr zeitaufwendig, da der Anwender die 
ganze Bildserie analysieren muss. Zum andern darf beim Aus- 
schneiden und Einfugen des segmentierten Organs im Ubergangs- 
bereich (Randbereich des Organs) keine Bildinf ormation verlo- 
ren gehen, was bei derzeitigen Segmentier-Algorithmen - die 
10 zudem (was Entwicklung und Rechenleistung betrifft) ausge- 
sprochen aufwendig sind - nicht garantiert wird. 

Aufgabe der vorliegenden Erfindung ist es daher, eine Technik 
zur Verbesserung einer organspezif ischen Bild-Optimierung 
15 vorzuschlagen. 

Erf indungsgemafi wird also ein Verfahren zur organspezif ischen 
Bildoptimierung in der Computer-Tomographie vorgeschlagen 
welches folgende Schritte aufweist: 

20 

- Aufnehmen einer Korperschicht mit einem CT-Gerat, 

- Berechnen der HU-Werte (Grau-Werte) dieser Schicht, 

- Darstellen der HU-Werte-Hauf igkeitsverteilung dieser 
Schicht in einem Histogramm, 

5 - Unterteilen des Histogramms in verschiedene (organspezif i- 
sche) HU-Bereiche, 

- Zuordnen von Bildfiltern unterschiedlicher Filterbreite zu 
den verschiedenen HU-Bereichen, 

- Filtern der zu dem jeweiligen HU-Bereich zugehorigen Bild- 
30 pixel mit dem entsprechendem Bildfilter. 

Dabei erfolgt das Zuordnen erf indungsgemaS durch eine HU- 
abhangige Transf erf unktion . 

3 5 In einer vorteilhaf ten Ausfuhrung der Erfindung wird zur Dar- 
stellung des Histogrammes ein vorgef iltertes Bild herangezo- 
gen welches zuvor aus dem CT-Bild erzeugt wurde. 



200201317 

» I 
• j 

7 

Dieses vorgef ilterte Bild ist entweder durch einen scharfen 
Bildkontrast oder durch ein geringes Bildrauschen ausgezeich- 
net . 

5 

Genauso kann es aber auch vorteilhaft sein aus dem erstellten 
CT-Bild zwei vorgef ilterte Bilder zu erzeugen, welche zur 
Darstellung des Histograinmes herangezogen werden, wobei eines 
der Bilder durch scharfen Kontrast und das andere Bild durch 
10 geringes Bildrauschen gekennzeichnet ist. 

0 Dabei werden erf indungsgemaS das oder die vorgef ilterten Bil- 
der bei der Bildrekonstruktion entweder durch einen entspre- 
chenden Faltungskern, durch ein zweidimensionales separierba- 
15 res Bildfilter oder durch jeweils zwei eindimensionale Bild- 
filter erzeugt. 

Bei der Verwendung zweier unterschiedlich vorgef ilterter Bil- 
der wird erf indungsgemaS ein optimiertes CT-Bild dadurch er- 

2 0 zeugt, dass beide vorgef ilterten Bilder mittels einer HU- 

abhangigen Transf erf unktion lokal gemischt werden. 

In einer weiteren Ausfiihrung der Erfindung werden die unter- 
schiedlich gefilterten HU-Bereiche vorteilhaf ter weise in un- 
*^ls terschiedlicher Fensterung dargestellt. 

Des weiteren wird gemafi der vorliegenden Erfindung ein Compu- 
ter-Tomographiegerat beansprucht, welches zur Durchfiihrung 
des Verfahrens nach einem der Anspriiche 1 bis 10 geeignet ist 

3 0 und einen Computer aufweist, der die einzelnen Schritte des 

erf indungsgemaSen Verfahrens durchfiihrt. 

Ferner wird ein Computersof tware-Produkt beansprucht, welches 
ein Verfahren nach einem der Anspriiche 1 bis 10 implemen- 
3 5 tiert, wenn es auf einer mit einem Computer-Tomographiegerat 
verbundenen Recheneinrichtung lauft. 
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Weitere Merkmale, Eigenschaf ten und Vorteile der vorliegenden 
Erfindung werden nunmehr anhand von Ausf iihrungsbeispielen und 
unter Bezugnahme auf die begleitenden Figuren der Zeichnungen 
erlautert . 

5 

Figur 1 zeigt schematisch eine CT-Apparatur fiir ein Facher- 
strahlverf ahren gemaS der vorliegenden Erfindung, 



Figur 2 zeigt eine Hounsf ield-Skala in der die Hounsfield- 
10 Einheiten (HE) unterschiedlicher Organe des menschli- 

chen Korpers angegeben sind, 

Figur 3 zeigt eine Fensterung (Windowing) bei der Darstellung 
von CT-Bildern, 

15 

Figur 4 zeigt ein Histogramm der HU-Werte im Thoraxbereich, 



Figur 5 zeigt eine Zuordnung unterschiedlicher diskreter Fil- 
terfunktionen (Transf erf unktionen) zu unterschiedli- 
2 0 chem anatomischen Gewebe, 



Figur 6 zeigt den Verlauf einer kontinuierlichen Transfer- 
funktion durch die CT-Bilder weicher und scharfer 
Bildcharakteristik (I S harp und Ismooth) gemischt werden, 

Figur 7 zeigt wie ein zweidimensionaler separierbarer Bild- 
filter durch zwei eindimensionale Bildfilter ersetzt 
werden kann . 



30 In Figur 1 ist schematisch ein Computer-Tomographiegerat fur 
ein Facherstrahlverf ahren dargestellt das gema£ der vorlie- 
genden Erfindung arbeitet. Bei diesem Gerat rotieren Rontgen- 
Rohre 1 und Strahlenempf anger 2 (Detektoren) gemeinsam um ei- 
ne Drehmitte, die auch Mitte des kreisf ormigen Messfeldes 5 

3 5 ist, und in der sich der zu untersuchende Patient 3 auf einer 
Patientenliege 4 befindet. Um verschiedene parallele Ebenen 
des Patienten 3 untersuchen zu konnen, kann die Patientenlie- 
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ge entlang der Korperlangsachse verschoben werden. Wie man 
aus der Zeichnung erkennen kann, ergeben sich bei CT- 
Aufnahmen Transversalschnittbilder , also Abbildungen von K6r~ 
perschichten, die im wesent lichen senkrecht zur Korperachse 
5 orientiert sind. Diese Schichtdarstellungsmethode stellt die 
Verteilung des Schwachungswertes m(*/Y) selbst dar (z ist 
die Position auf der Korperlangsachse) . Die Computer- 
Tomographie (im folgenden CT genannt) benotigt Projektionen 
unter sehr vielen Winkeln a. Zur Erzeugung einer Schichtauf- 

10 nahme wird der von der Rontgenrohre 1 emittierte Strahlenke- 
gel so ausgeblendet , dass ein ebener Strahlenf acher entsteht, 
^ der eindimensionale Zentralpro j ektionen der durchstrahlten 
Schicht entwirft. Zur exakten Rekonstruktion der Verteilung 
der Schwachungswerte li z (x,y) muss dieser Strahlenf acher senk- 

15 recht auf der Drehachse stehen und aufierdem so weit gespreizt 
sein, dass er aus jeder Projektionsrichtung a die anvisierte 
Schicht des Messobjektes vollstandig iiberdeckt . Dieser das 
Objekt durchdringende Strahlenf acher wird von Detektoren die 
auf einem Kreissegment linear angeordnet sind aufgefangen. 

20 Bei handel sub lichen Geraten sind dies bis zu 1000 Detektoren. 
Der einzelne Detektor reagiert auf die eintref f enden Strahlen 
mit elektrischen Signalen, deren Amplitude proportional zur 
Intensitat dieser Strahlen ist. 

Jedes einzelne zu einer Projektion cc gehorige Detektorsignal 
wird jeweils von einer Messelektronik 7 aufgenommen und an 
einen Computer 8 weitergeleitet . Mit dem Computer 8 lassen 
sich die gemessenen Daten nun in geeigneter Weise verarbeiten 
und zunachst in Form eines Sinugramms (in dem die Projektion 
30 a als Funktion der Messwerte des entsprechenden Kanals P auf- 
getragen wird) in sogenannten Gordon-Einheiten, schlieiSlich 
aber in Form eines natiir lichen Rontgenbildes in Houns field- 
Einheiten an einem Monitor 6 visualisieren . 

35 Ziel der vorliegenden Erfindung ist es, durch ein geeignetes 
Verfahren, das letztlich im Computer 8 implementiert und 
durchgefuhrt werden soil, auf einfache und schnelle Weise ein 
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CT-Bild zu konstruieren in dem die unterschiedlichen anatomi- 
schen Gewebstypen mit unterschiedlicher Bildcharakteristik 
(Bildscharfe und Bildrauschen) dargestellt sind und somit ein 
organspezif isch optimiertes CT-Bild erhalten wird. 

5 

Die vorliegende Erfindung macht sich die Tatsache zunutze, 
dass die CT-Werte (HU-Werte) je nach Organ in der Hounsfield- 
Skala unterschiedliche Bereiche belegen. 

10 Die Hounsf ield-Skala ist in Figur 2 dargestellt. Die CT-Werte 
der einzelnen Organe nehmen weitgehend unabhangig vom verwen- 
p deten Rontgenspektrum spezifische Bereiche ein. So weisen 
Lungengewebe und Fett wegen ihrer niedrigen Dichte und der 
dadurch bedingten niedrigen Schwachung negative CT-Werte im 

15 Bereich von -950 bis -550 bzw. -100 bis -80 auf . Die meisten 
anderen Gewebe liegen im positiven Bereich (Niere: 20-40, 
Herz: 40-100, Blut: 50-60, Leber: 50-70). Knochengewebe be- 
sitzt wegen der hohen Ordnungszahl von Kalzium und damit der 
hoheren Schwachung CT-Werte bis zu 2000 HU. 

20 

Fur herkommliche CT-Gerate stehen zur Darstellung der gesam- 
ten Hounsf ield-Skala 4096 (=2 12 ) unterschiedliche Grauwerte 
zur Verfugung. Vom Betrachter konnen jedoch nur maximal 60 
bis 80 Graustufen unterschieden werden. Aus diesem Grund wird 
bei der CT-Bildgebung - wie bereits erwahnt - eine Fensterung 
(engl.: Windowing) vorgenommen in dem einem interessierenden 
HU-Intervall die gesamte Grauwerteskala zugeordnet wird. In 
Figur 3 ist dies beispielsweise fur den Bereich kompakter 
Knochensubstanz dargestellt. Interaktiv wird vom Anwender das 
30 Fenster durch Zentrum (Fensterlage C) und Weite (Fensterbrei- 
te W) z.B. per Maus oder Drehknopf festgelegt. Im Beispiel 
der Figur 3 liegt das Zentrum bei C=2 000, die Fensterbreite 
bei W=400. Dabei sind dem Fenster zwischen weifi und schwarz 
10 Graustufen zugeteilt. 



35 



Der erf inderische Gedanke besteht nun darin, ein CT-Bild or- 
ganspezif isch auf Basis der CT-Werte in der HU-Skala und 
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durch entsprechende Fensterung zu optimieren. Erf indungsgemalS 
wird vorgeschlagen, zunachst die CT-Werte (HU-Werte) des Bil- 
des einer aufgenommenen Schicht in einem Histogramm darzu- 
stellen. Ein solches Histogramm ist beispielhaft fur den Tho- 
5 raxbereich in Figur 4 dargestellt. Einzelne HU-Bereiche des 
Histogramms (A, B, C, D, E, F) korrelieren mit spezifischen Organ- 
strukturen - oder anders ausgedruckt: Segment-Abschni tte 
(A, B, C, D, E, F) der HU-Verteilungskurve 9 sind eindeutig einem 
Organ bzw. einer Gewebsstruktur zugeordnet. Im Falle des 
10 Histogrammes der Figur 4 korreliert der HU-Bereich A bei- 

spielsweise mit Lungengewebe , der Bereich B mit Fettgewebe, C 
mit eingelagertem Wasser, D mit Blut (z.B. Aorta), E mit Le- 
bergewebe und F mit Herzgewebe. 

15 Urn eine organspezif ische Einstellung der Bildcharakteristik 
(Bildrauschen und Bildscharfe) vornehmen zu konnen wird im 
einfachsten Fall der erf indungsgemafien Ausflihrung jedem or- 
ganspezif ischen Kurvenabschnitt HU(x,y) eine eigene Transfer- 
funktion A,(HU(x,y)) zugeordnet, X kann als Filterf unktion 

20 aufgefasst werden mit der samtliche Bildpixel I(x,y) des je- 
weiligen HU-Bereiches gefiltert werden. Mathematisch lasst 
sich dies folgendermaSen darstellen: 

HMHU(x 9 y)) = MHU(x 9 y))-I(x,y) 

Im einfachsten Fall ist X eine konstante Funktion die fur je- 
den organspezif ischen HU-Bereich einen speziellen Wert hat. 
Dieser spezielle Wert ist vom Anwender unter Berucksichtigung 
der erwlinschten Bildcharakteristik des entsprechenden Organs 
30 festzulegen. Dies fuhrt dazu, dass im CT-Bild die Pixelberei- 
che die zu entsprechenden organspezif ischen Grauwertzonen ge- 
horen unter schiedlich gefiltert bzw. manipuliert werden. 

Figur 5 stellt vereinfacht ein aus den Rohdaten rekonstruier- 
35 tes CT-Bild 10 dar, in dem die organspezif ischen Grauwertzo- 
nen, die den Histogramm-Abschnitten aus Figur 4 zugehoren, 
eingekreist sind. Wie man am HU-Bereich A sieht, mussen die 
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Grauwertzonen nicht unbedingt zusammenhangen. Im Beispiel der 
Figuren 4 und 5 werden die organspezif ischen Bereiche 
A,B,C,D,E,F HU-basiert gefiltert. 

Diese HU-basierte Filterung lasst sich mathematisch folgen- 
dermafien beschreiben: 



7(4) = ^-IiHUix.y)), HUeA 
/(^) = A.IiHUix.y)), HUeB 

IW = \-I(HU(x,y)), HUeF 



Erf indungsgemaS wird jede Transf erf unktion Xi so gewahlt, 
dass das entsprechende Organ bzw. Gewebe i in der gewunschten 
15 optimalen Bildcharakteristik erscheint. 

Nach einer Bildf ilterung gemaS obiger Vorgehensweise werden 
in einem weiteren Schritt der Erfindung die unterschiedlichen 
HU-Bereiche in unterschiedlicher Fensterung (Windowing) dar- 
2 0 gestellt. Dabei wird als Zentrum des Fensters der mittlere 

HU-Wert der interessierenden Struktur i gewahlt. Die Fenster- 
breite richtet sich nach den Schwachungsdif ferenzen der je- 
weiligen Struktur: Fur die Darstellung sehr kleiner Schwa- 
chungsunterschiede, wie z.B. bei Hirngewebe, wird ein enges 
5 Fenster, bei grolSen Schwachungsunterschieden, wie z.B. bei 
Lungengewebe , wird ein weites Fenster gewahlt. 



Der Vorteil des erf indungsgemaSen Verfahrens ist anhand Figur 
5 deutlich zu sehen: Nach dem Verfahren im Stand der Technik 

3 0 miissen hinsichtlich der Bildcharakteristik eine Vielzahl von 
CT-Bildern 10 erzeugt, von Hand selektiert und in den ausge- 
wahlten CT-Bildern dann die entsprechenden Organs truktur en 
segment iert und in das endgiiltige Bild eingeftigt werden. In 
der vorliegenden Erfindung vermittelt das automatisch er- 

3 5 stellte Histogramm der Bilddaten die Beziehung zwischen der 
Organstruktur i und den zugeordneten Bildelementen (Pixeln: 
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x,y) . Es muss lediglich die jeweilige Transf erf unktion X± vom 
Anwender vorgegeben werden. 

Unter Umstanden ist es wunschenswert das Histogramm erst nach 
5 einer zuvor durchgef iihrten Manipulation des rekonstruierten 
CT-Bildes in Form einer Bildf ilterung zu erstellen. Bei- 
spielsweise gelingt es nach einer Glattung des CT-Bildes Ge- 
faSe im Lungenparenchym mit der gleichen Scharfe wie das Lun- 
genparenchym selbst darzustellen (Ein Parenchym ist der Zel- 
10 lenverband eines Organs der dessen Funktion bestimmt) . Wird 

das CT-Bild einer Lungenaufnahme hinreichend stark geglattet, 
so kann die gesamte Lunge einschliefilich deren (Lungen- 
fremden) Gefafistrukturen in einem zusammenhangenden Intervall 
des Histogrammes abgebildet werden. Ein derart geglattetes 
15 Ausgangsbild wird im folgenden als I W eich (engl.: I S mooth ) be- 
zeichnet. Genauso aber kann es sein, dass zur Histogram- 
Erstellung ein kontrastreiches Ausgangsbild I S charf (engl.: I s _ 
harp) erforderlich ist um beispielsweise auEerorganisches Ge- 
webe mit ebenso hoher Scharfe darzustellen. 

20 

Die jeweilige HU-basierte Bildf ilterung lasst sich in diesen 
beiden Fallen f olgendermaSen darstellen: 

7(4) = A,-I weich (HU(x,y)) 9 HUei 

5 bzw. 

M) = X r I scha j(HU{x,y)), HUei 



In einer weiteren Ausfiihrung der Erfindung ist die Transfer- 
funktion X, mit der die organspezif ischen Histogrammbereiche 

30 gefiltert werden, nicht diskret, sondern hat - wie beispiels- 
weise in Figur 6 dargestellt - einen kontinuierlichen Verlauf 
iiber den gesamten Histogrammbereich. Der Verlauf wird vom An- 
wender fur jede zu messende Korperschicht festgelegt. Der Be- 
trag von X bewegt sich in einem Intervall zwischen 0 und 1. 

3 5 Eine derartige HU-basierte kontinuierliche Bildf ilterung 

lasst sich mathematisch als eine Mischung zweier extremaler 
CT-Ausgangsbilder I W eich und I SC harf beschreiben: 
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I(A(HU(x 9 y)) = MHU(x,y))-I schaff (x,y) + (l-MHU(x,y)))I weich (x 9 y) 

Die Erzeugung von CT-Ausgangsbildern I SC harf bzw. Iweich / also 
5 von CT-Bildern mit extrem scharfer bzw. extrem weicher Bild- 
charakteristik, lasst sich auf unterschiedliche Weise reali- 
sieren. Bereits bei der Rekonstruktion des Bildes aus den CT- 
Rohdaten - beispielsweise durch Fouriertransf ormation - ist 
es moglich die Bildcharakteristik vor der Rucktransf ormation 

10 durch Wahl eines entsprechenden Faltungskerns (bzw. dessen 

Fouriertransf ormierte) zu beeinf lussen . Ein scharfer kanten- 

^ betonender Faltungskern liefert ein kontrastreiches Ausgangs- 
bild Ischarf bei entsprechend hohem Bildrauschen wahrend ein 
weicher glattender Faltungskern zwar die Ortsauf losung redu- 

15 ziert jedoch ein rauscharmes Bild I W eich erzeugt. 

Eine andere Moglichkeit besteht darin, die durch eine CT- 
Faltungskern hervorgeruf ene Bildcharakteristik durch einen 
zweidimensionalen separierbaren Bildfilter retrospektiv aus 
20 dem bereits rekonstruierten CT-Ausgangsbild einzustellen. 

Noch einfacher ist es einen zweidimensionalen separierbaren 
Bildfilter durch zwei eindimensionale Bildfilter zu ersetzen. 

Wie in Figur 7 dargestellt ist werden bei einem zweidimensio- 
5 nalen separierbaren Bildfilter mit einer 1-Pixelbreite zur 
Filterung eines Pixels insgesamt 9 Pixel miteinbezogen . Wird 
der zweidimensionale separierbare Bildfilter gemaS Figur 7 
durch zwei eindimensionale Bildfilter mit jeweils einer 1- 
Pixelbreite ersetzt, so werden bei der Filterung pro Pixel in 
3 0 jeder Dimension nur 3 Pixel, insgesamt also 6 Pixel, mitein- 
bezogen . 

Letztendlich bedeutet eine retrospektive Manipulation des CT- 
Bildes mit zwei eindimensionalen Bildfiltern zum einen eine 
35 deutliche Reduzierung der Rechenzeit. Zum andern macht ein 

Filtern des rekonstruierten Bildes das Speichern der Rohdaten 
iiberf liissig. 
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Das Verfahren bzw. die Technik um die durch einen CT- 
Faltungskern hervorgeruf ene Bildcharakteristik durch einen 
zweidimensionalen separierbaren bzw. durch zwei eindimensio- 
nale Bildfilter einzustellen ist in T.Flohr, S.Schaller, 
A.Stadtler, et al . , "Fast image filters as an alternative to 
reconstruction kernels in Computed Tomography", Proceedings 
of SPIE, Vol.4322 (2001), pp 924-933 ausfiihrlich dargestellt. 
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Patentanspriiche 

1. Verfahren zur organspezif ischen Bildoptimierung in der 
Computer-Tomographie aufweisend die folgenden Schritte : 

5 

- Aufnehmen einer Korperschicht mit einem CT-Gerat, 

- Erstellen eines CT-Bildes dieser Schicht durch Berechnen 
der HU-Werte (Grau-Werte) dieser Schicht, 

- Erstellen einer HU-Werte-Hauf igkeitsverteilung des CT- 
10 Bildes dieser Schicht in einem Histogramm, 

- Unterteilen des Histogramms in verschiedene organspezif i- 
^ sche HU-Bereiche, 

- Zuordnen von Bildfiltern unterschiedlicher Filterbreite zu 
den verschiedenen HU-Bereichen, 

15 - Filtern der zu dem jeweiligen HU-Bereich zugehorigen Bild- 
pixel mit dem entsprechendem Bildfilter. 

2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, dass das 

20 Zuordnen durch eine HU-abhangige Transf erfunktion erfolgt. 

3. Verfahren nach Anspruch 1 bis 2, 

dadurch gekennzeichnet, dass aus 

dem erstellten CT-Bild ein vorgef iltertes Bild erzeugt wird 
-*%5 welches zur Darstellung des Histogrammes herangezogen wird. 

4. Verfahren nach Anspruch 1 bis 2, 

dadurch gekennzeichnet, dass das 

vorgef ilterte Bild durch einen scharfen Bildkontrast ausge- 
30 zeichnet ist. 

5. Verfahren nach Anspruch 1 bis 2, 

dadurch gekennzeichnet, dass das 

vorgef ilterte Bild durch ein geringes Bildrauschen ausge- 
35 zeichnet ist. 
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6. Verfahren nach Anspruch 1 bis 3, 

dadurch gekennzeichnet, dass aus 

dem erstellten CT-Bild zwei vorgef ilterte Bilder erzeugt wer- 
den welche zur Darstellung des Histogrammes herangezogen wer- 
5 den, wobei eines der Bilder durch scharfen Kontrast und das 
andere Bild durch geringes Bildrauschen gekennzeichnet ist. 



7. Verfahren nach Anspruch 1 bis 6, 

dadurch gekennzeichnet, dass das 

10 oder die vorgef ilterten Bilder bei der Bildrekonstruktion 
durch einen entsprechenden Faltungskern erzeugt werden. 

► 

8. Verfahren nach Anspruch 1 bis 6, 

dadurch gekennzeichnet, dass das 

15 oder die vorgef ilterten Bilder bei der Bildrekonstruktion 

durch ein zweidimensionales separierbares Bildfilter erzeugt 
werden . 



9. Verfahren nach Anspruch 1 bis 6, 

20 dadurch gekennzeichnet, dass das 

oder die vorgef ilterten Bilder bei der Bildrekonstruktion 
durch jeweils zwei eindimensionale Bildfilter erzeugt werden. 

10. Verfahren nach Anspruch 1 bis 6, 

Jfk5 dadurch gekennzeichnet, dass aus 

den zwei unterschiedlich vorgef ilterten Bildern ein optimier- 
tes CT-Bild erzeugt wird, indem beide vorgef ilterten Bilder 
mittels einer HU-abhangigen Trans f erf unkt ion lokal gemischt 
werden . 

30 

11. Verfahren nach Anspruch 1 bis 10, 

dadurch gekennzeichnet, dass die 

unterschiedlich gefilterten HU-Bereiche in unterschiedlicher 
Fensterung dargestellt werden. 



35 
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12 . Computer-Tomographiegerat zur Durchfuhrung des Verf ahrens 
nach einem der Anspriiche 1 bis 10, mit einem Computer, der 
die einzelnen Schritte der Signalverarbeitung durchfiihrt. 

13 . Computersof tware-Produkt , 

dadurch gekennzeichnet, dass es 

ein Verfahren nach einem der Anspriiche 1 bis 10 implemen- 
tiert, wenn es auf einer mit einem Computer-Tomographiegerat 
verbundenen Recheneinrichtung lauft. 
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Zusammenf assung 

Verfahren zur Grauwert-basierten Bildf ilterung in der Compu- 
ter -Tomographie 

5 

Die vorliegende Erfindung bezieht sich auf ein Verfahren zur 
organspezif ischen Bildoptimierung in der Computer-Tomographie 
aufweisend die folgenden Schritte : 

10 - Aufnehmen einer Korperschicht mit einem CT-Gerat, 

- Erstellen eines CT-Bildes dieser Schicht durch Berechnen 
der HU-Werte (Grau-Werte) dieser Schicht, 

- Erstellen einer HU-Werte-Hauf igkeitsverteilung des CT- 
Bildes dieser Schicht in einem Histogramm, 

15 - Unterteilen des Histogramms in verschiedene organspezif i- 
sche HU-Bereiche, 

- Zuordnen von Bildfiltern unterschiedlicher Filterbreite zu 
den verschiedenen HU-Bereichen, 

- Filtern der zu dem jeweiligen HU-Bereich zugehorigen Bild- 
20 pixel mit dem entsprechendem Bildf ilter. 

Fig. 4 
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